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Abstract 

A phenomenological theory of electric double layer polarization of blocking electrode is presented solving modified 
Debye-Falkenhagen (MDF) equation for potential, under the impedance boundary condition. The dynamic impedance 
and capacitance are obtained in terms Debye screening length and frequency dependent polarization length. Two char- 
acteristic relaxation frequencies for compact layer and diffuse layer are identified. At frequencies less than Helmholtz 
layer relaxation frequency the EDL is perfectly blocking. At frequencies larger than the diffuse layer relaxation fre- 
quency the EDL behaves like a resistor. At crossover frequencies the electrode is polarized. Theoretical results agrees 
well with experimental capacitance dispersion data. 
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1. Introduction 

Electrode polarization |fl S ! & 0, & is a ubiq- 
uitous phenomenon and is found to have an increasing ap- 
plication in supercapacitors ifioll . batteries ifTHl . fuel cells 
11211 . solar cells Ml 311 and ion desalination membrane elec- 
trodes [14]. Relaxation and dispersion lfl5l [l6 , 17 1 are 
two electrochemical phenomena usually observed when 
the electrode is coupled to various electrolytes viz., aque- 
ous, liquid crystals glassy electrolytes ifioh - polyelec- 
troly tes lfl8ll and ionic liquids fl5ll . 

In this communication, we report a model of electrode 
polarization in EDL formed near a planar blocking elec- 
trode (Fig.l). The EDL is divided into two region, com- 
pact layer and diffuse layer separated by outer Helmholtz 
plane (OHP). The objective of the paper is to understand 
the interplay between the compact and diffuse layers in- 
terns of the length scale and various time scales emerging 
out of the electric double layer phenomena when a sinu- 
soidal time varying potential is applied. 



2. Model and Solution 

The Debye-Falkenhagen equation |6, 8, 19] for relax- 
ation of the ionic atmosphere in diffuse layer at the planar 
electrode is 

1 hn . 

(1) 



D dt H 



where D is the diffusion coefficient, p - e{z + c + - Z-cJ) 
is the charge density and c ± is th ionic concentrations, z± 
the ionic charges and k~ x is the Debye screening length 
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where / = (1/2) £" =1 c,z 2 is the ionic strength, k% the 
Boltzmann's constant, Na is tha Avogadro number, T the 
temperature, eo the permittivity of free space and e r di- 
electric constant of the solvent. Now we assume near 
equilibrium local concentration ions c,- is given by Boltz- 
mann equation 



d = c°° exp(±Zie4>/k B T) 



(3) 
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Now for a symmetrical z+ — z~ — z electrolytes and as- 
suming c+ = c!° = co, linearizing Eq. (0 under condition 
<p «: ksT/ze and substituting in Eq. (|T), we obtained the 
modified Debye-Falkenhagen (MDF) equation for poten- 
tial d> 



D dt Y 



(4) 



For a sinusoidal applied potential, cb - <f>(f)e l0)t the MDF 
equation is 



[V 2 -k 2 (1 +iiOT D )] 0(r) = 



(5) 



where tq = 1/k 2 D is the Debye time 12011 . Eq. Q rep- 
resent the dynamic picture of the interface, when a sinu- 
soidal potential is applied at the interface (Fig. 1). 
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Figure 1 : Schematic diagram of electric double layer (EDL) formed at 
planar blocking electrode. 

In the case of blocking electrode, the ohmic current (J) 
crossing a metal/electrolyte interface is given by 



-crV0 



(6) 



where cr is specific interfacial conductivity of electrolyte 
near electrode. The interfacial and bulk electrolyte con- 
ductivity ( o"o) satisfies cr » cr lUsh . where <x = e e r /cD 
lHHD]. Eq. (|6]l gives the current arriving at the OHP and is 
equal to the ratio of potential drop at the OHP (<p - c}), to 
the impedance (Z#) of the Helmholtz layer, i.e. (</>-0)/Z# 
where (f> is the potential of the electrode and <p is potential 
at the OHP. Now equating the current on both sides we 
have the impedance boundary condition 02 ill 



-cr(d<l>/dz) = (<p-$)/Z H 



(7) 



where Z# = (l/icoCn + Rh) is the impedance of the 
Helmholtz layer and Rh and Ch is the compact layer re- 
sistance and capacitance. The other semi-infinite (bulk) 
boundary condition is (f>(oo) = 0. The impedance 
boundary condition may also be written as <p = - 
hc{(i>y x (d(j) I dz) where A c (w) = crZ H . The quantity 
| Ac(<w) r 1 is the phenomenological frequency depen- 
dent length. It is the effective thickness of interfacial layer 
where the electrode is polarized and decides the crossover 
time scales from purely capacitive regimes to resistive 
regimes we called it the "polarization length". It may 
be rewritten for simplicity as Ac(aS) = l/(o-/ia>CH + W) 
where W = crR H is a characteristic resistive length which 
may be called as modified Wagner number 112.21 l23ll . The 
value of [ Ac(w) | _1 is decided by the compact layer ca- 
pacitance Ch, resistance Rh and the conductivity of the 
electrolyte cr. 

We use the method developed in ref. ll24ll to obtained 
solution of MDF Eq. ©. The total admittance of planar 
blocking electrode is 



Yp(ai) = 



Z H + l/(crK(l+iajT D y' 2 ) 



(8) 



where An, is the area of electrode. In order to understand 
behavior of admittance, we consider the following 
frequency regimes. 

At low frequency /e~' <k| Ac(oj) from Eq. ([8]l we 
have 



j _ AcM | / A c (fa/n _^ 



(9) 



The leading contribution in admittance originates from 
the Helmholtz layer and its characteristic relaxation 
frequency is u>h = (ChRh) • 1 ■ When the frequency is 
co <K u>h, Eq. © shows that the EDL behaves like an 
ideally polarizable blocking electrode. 

At high frequency | Ac(w) T 1 ^ K~ l , from Eq. (IS) we 
have 



Yp(cS) = Aqctk 



1 - 
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When the frequency is higher than the characteristic dif- 
fuse layer relaxation frequency, cap, < u, Eq. (TlOb shows 
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that the leading order contribution to the interfacial ad- 
mittance comes from the diffuse layer resistance. This 
observation will allow us to see the effects of concentra- 
tion on the admittance of planar electrode. The admit- 
tance is resistive and the characteristic frequency at which 
the EDL relaxes is the characteristic frequency at which 
diffuse layer relaxes, too = (ctk)/Cd = <t/(6o€d). But 
ct/(£q€o), is also a characteristic conductivity relaxation 
frequency [ 1 ] . 

At the crossover region, the frequency range is confined 
to loh < to < lod and | Ac(to) \ is comparable to k . 
In this case the crossover frequency is resultant of both 
processes- relaxation of Helmholtz layer and the relax- 
ation of diffuse layer; and we expect the dispersions in 
electrochemical response of EDL. 

Now expanding Eq. ([8]) in the ratio of the Debye screen- 
ing length to ion diffusion length, we have 

1 i I to \ 2 

Y P (to) = Aqklt H — — 

[1+kltZh 2(1 +KtrZ H ) 2 \k 2 D) 

1 + 3kctZ h I to \ 4 

+ %{1+kltZ h ? \k 2 DI 

At high frequencies, k~ 1 / y/D/to > 1, then the higher 
terms will contribute. Gunning et al [7] obtained the 
impedance depended on zeta potential of the diffuse dou- 
ble layer. The theory developed here differs in the inclu- 
sion of the Stern layer along with the diffuse double layer. 
Our theory focuses on the boundary constraint at the OHP 
rather than the "slipping plane" H. 

3. Results and Discussions 

Fig. 2 (a) show the effect of concentration on the log- 
log plot of magnitude of impedance, \Z(to)\ and frequency 
to. As we increase the concentration of electrolyte the 
crossover frequency (system changes from capacitive to 
resistive behavior) decreases. The frequency dependence 
of electrode polarization is shifted to lower frequencies as 
the concentration of the electrolyte decrease. Fig. 2 (b) 
shows the phase angle <t> vs log to. The plot shows that 
the impedance behavior is resistive at the high frequency, 
which is identified in the plot as cD — » 0°. The capacitive 
behavior is seen at a low frequency which is identified 
in the plot as O — > 90°. As we decrease the concentra- 
tion of electrolyte the crossover frequency at which the 
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Figure 2: Effect of concentration variation is shown in (a) and (b). The 
Debye length k~ 1 is varied from 0.3, 1 and 3. Effect of electrolyte con- 
ductivity variation is shown in (c) and (d). The diffusion coefficient 
D (10~ 5 cm 1 Is) is varied from 0.02, 1 and 4. The graphs are plotted 
using Eo = 8.854 xl0~ 14 (F/cm), en = 3.5, ro= 0.53 nm and area of 
electrode A = 0.038cm 2 



phase behavior switch from capacitive to resistive behav- 
ior is shifted progressively to lower frequencies. Fig. 2(c) 
show the influence of electrolyte conductivity through the 
diffusion coefficient of ions on the log-log plot of mag- 
nitude of impedance, \Z(co)\ and frequency to. With the 
increase in the diffusion coefficient of electrolyte the fre- 
quency at which the system changes from capacitive to 
resistive behavior increases. The frequency dependence 
of electrode polarization occurs at a higher frequency as 
we increase the conductivity (by increasing the diffusion 
coefficient). Thus in viscous solutions where the diffusion 
coefficient of ion is low the frequency dependent behav- 
ior is shifted progressively to lower frequencies. Fig. 2 
(d) shows the phase angle O vs log to. The plot shows that 
the impedance behavior is resistive at the high frequency, 
which is identified in the plot as O — > 0°. The capacitive 
behavior is seen at a low frequency which is identified 
in the plot as O — > 90°. As we increase conductivity of 
electrolyte the characteristic frequency at which the phase 
behavior crossover from capacitive to resistive behavior is 
shifted progressively to higher frequencies. 
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Fig. 3 shows the theoretical plots fits with experimen- 
tal capacitance dispersion data for Hg/electrolyte inter- 
face. The real part of the capacitance is calculated as 
C'(a>) - Z" (oj) / aj\Z(to)\ 2 . In agreement with the theory, 
the data shows two crossover frequencies. The effect of 
electrolyte conductivity on capacitance is exactly same 
as in experiment. Increasing the electrolyte conductivity 
by increasing diffusion coefficient, increases the crossover 
frequency {cdd) of capacitance dispersion. 
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Figure 3: Comparison of theory with experiment. Data were taken from 
ref. tl6Lll7ll . Physical parameters used to fit the data are listed in table 
1. 



4. Conclusions 

The electrode polarization is strongly dependent on the 
phenomenological polarization length. The Helmholtz 
and diffused layers relaxation frequencies divide the 
whole frequency range into three regimes. The low fre- 
quency regime is characterized by large polarization layer 
thickness and the capacitance of EDL is that of Helmholtz 
layer. Whereas the high frequency is characterized by 
small polarization layer and the capacitance of the EDL is 
that of diffuse layer. The crossover frequency is indicated 
by a transition from capacitance to resistive behavior and 
an S-shaped crossover region in the capacitance plot. The 
capacitance dispersion that arises in the crossover regime 
is found to depend on the concentration and conductiv- 
ity of the electrolyte. Decrease in the concentration of 
electrolyte results in the shift of the crossover region pro- 
gressively to lower frequencies. The capacitance disper- 
sion is seen for frequencies where % < u> < o>d- The 



Table 1 : Here is e# dielectric constant and r# (nm) is effective thickness 
(corresponding to hydrated ion size) of Helmholtz layer. The effective 
diffusion coefficient D ( 10~ 6 cm 2 /sec) used as fitting values are listed. 
Here we have assumed that two conductivities follow: cr = kcrg 11511 

and k = 1CT 6 . 

Data E S C b e H £g rg A D 

Fig.5(a) NH 4 F H 2 1(T 3 3.5 78.6 0.53" 38 18.3 
Fig.5(b) LiCl gly 3.2 6 42.5* 0.37" 25 0.06 

Fig.5(b) LiCl gly+H 2 Q 3.2 6 42.5* 0.37" 25 15 
E- Electrolyte, S- solvent.gly- Glycerol, Cb- molar concentration, Aq 
(10~ 3 )- area of electrode in cm 2 , a- [25], b- 1 26] 

theory agrees with the capacitance dispersion data of the 
mercury/electrolyte interface for an aqueous, glycerol and 
glycerol-water system. Finally one can say, this theory is 
an indispensable step in the quantitative description of the 
dynamics of planar blocking electrode. 
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